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ABSTRACT 

Observations of the hot gas, which is surrounding Sgr A* and a few other nearby 
galactic nuclei, imply that mean free paths of electron and proton are comparable to 
gas capture radius. So, hot accretion flows likely proceed under week collision con- 
ditions. As a result thermal conduction by ions has a considerable contribution in 
transfer of the realized heat in accretion mechanisms. We study a 2D advective ac- 
cretion disk bathed in a poloidal magnetic field of a central accretor in the presence 
of thermal conduction. We find self-similar solutions for an axisymmetric, rotating, 
steady, viscose-resistive, magnetized accretion flow. The dominant mechanism of en- 
ergy dissipation is assumed to be turbulence viscosity and magnetic diffusivity due to 
magnetic field of the central accretor. We show that the global structure of ADAF s 
are sensitive to viscosity, advection and thermal conduction parameters. We discuss 
how radial flow, angular velocity and density of accretion flows may vary with the 
advection, thermal conduction and viscous parameters. 
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1 INTRODUCTION 

The foundations of our present understanding about advec- 
tion dominated accretion flows were laid out in a series of 
papers by Narayan & Yi (1994,1995 a,b), although some 
ideas were anticipated much earlier by Ichimaro (1977). The 
specific abbreviation ADAF , which stands for advection 
dominated accretion flow , was introduced by Lasota et al. 
(1996). An ADAF is defined as one in which a large frac- 
tion of the viscously generated heat is advected with the 
accreting gas, and only a small fraction of the energy is ra- 
diated. ADAFs have an opposite regime in comparison with 
the standard model. In the standard model, flow is described 
in away that the heat generated by the viscosityradiates out 
of the system immediately after its generation (Shakura & 
Sunyaev 1973). 

These advection-dominated accretion flows occur in 
two regimes depending on their mass accretion and optical 
depth. Actually the optical depth of accretion flows is highly 
dependent on accretion rates. In a high mass accretion rate, 
the optical depth becomes very high and the radiation gen- 
erated by the accretion flow can be trapped within the disk. 
This type of accretion disks is known as optically thick or 
Slim Disks which has been introduced by Abramowicz et al. 
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(1998). In the limit of low mass accretion rate, the disk be- 
comes optically thin. In this case, the cooling time of accre- 
tion flows is longer than accreting time scale. So the energy 
generated by accretion flows mostly remains in the disks and 
the disks can not radiate their energies efficiently. This kind 
of accretion flows are named radiation inefficient accretion 
flows (RIAFs). This type of accretion flows is investigated 
by many authors (Narayan & Yi 1994, Abramowicz et al. 
1995, Chen 1995). 

These types of solutions have been used to interpret 
the spectra of X-ray binary black holes in their quiescent 
or low/hard state as an alternative to the Shapiro, Light- 
man and Eardly (1976, SLE) solutions. Since ADAFs have 
large radial velocities and also the infalling matter caries the 
thermal energy to the black hole, the energy transported by 
advection can stabilize the thermal instability by removing 
their steep temperature gradient; thus the ADAF models 
have been widely used to explain the observations of low lu- 
minous observed in Sgr A* (Narayan et al. 1996, Hameuray 
et al. 1997). However, numerical simulations of radiation in- 
efficient accretion flows revealed that the low viscous flows 
are convectively unstable and therefore convection strongly 
influences the global structure of accretion flows (Igumen- 
shchev, Abramowicz & Narayan 2000). Thus, another type 
of accretion flows was proposed, in which the convection 
plays a dominat mechanism in transporting the energy, an- 
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gular momentum and local released viscose energy within 
the disk. 

A remarkable problem arises when the accretion disks 
threaded by a magnetic field. In the ADAF models, the tem- 
perature of accreting disks is so high that the accreting ma- 
terials are ionized. So, the magnetic field plays an important 
role in the dynamics of the accretion flows. Some authors 
tried to solve the MHD equations of magnetized ADAFs an- 
alytically. For example, Kaburaki (2000) has presented a set 
of analytical solutions for a fully advective accretion flow in 
a global magnetic field. Shadmehri (2004) has extended this 
analysis for a non-constant resistivity. Ghanbari et al. (2007) 
have presented a set of self-similar solutions for 2D viscous- 
resistive advection dominated accretion flows (ADAFs) in 
the presence of dipolar magnetic field of the central accre- 
tor. They have shown that the presence of magnetic field 
and its associated resistivity can considerably change the 
picture of the accretion flows. 

Recent observations of hot accretion flows around active 
galactic nuclei indicate that they should be based on the 
collisionless regimes. 

Chandra observations provide tight constraints on both 
density and temperature of gas at or near the Bondi cap- 
ture radius in Sgr A* and several other nearby galactic nu- 
clei. Tanaka & Menou (2006) with some calculations have 
shown that the accretion disks in such systems will proceed 
under weakly-collisional conditions. So, the thermal conduc- 
tion has an important role in the energy transport along 
the disks. The aim of this work is to consider the effect of 
thermal conduction which has been largely neglected before, 
as an energy transport mechanism, on the 2D structure of 
ADAFs. It could affect the global properties of hot accre- 
tion flows substantially. A few authors considered the role 
of turbulent heat transport in ADAF disks (Honma 1996 , 
Manmoto et al.2000). Since thermal conduction acts to op- 
pose the formation of temperature gradient that causes it, 
one might expect that the temperature and density profiles 
for accretion flows are modified in which thermal conduc- 
tion plays a significant role to appear different, compared to 
those flows which thermal conduction is less effective (Shad- 
mehri 2008). 

The weakly-collisions nature of hot accretion flows has 
been addressed previously (Mahadevan & Quataret 1997). 
Johanson & Quataret (2007) studied the effect of elec- 
tron thermal conduction on the properties of hot accretion 
flows under the assumption of spherical symmetry. In an- 
other interesting analysis, Tanaka & Menou (2006), stud- 
ied the effect of saturated thermal conduction on optically 
thin ADAFs using an extension of self-similar solution of 
Narayan & Yi (1994). In their solutions, the thermal con- 
duction is provided an extra degree of freedom which affects 
the global dynamical behaviors of the accretion flow. Ab- 
bassi et al. (2008) have presented a set of self-similar solu- 
tions for ADAFs with a toroidal magnetic field in which the 
saturated thermal conduction has a great role in the energy 
transport in the radial direction. The tangled magnetic field 
in accretion flows would likely reduce the effective mean free 
paths of particles. The magnitude of this reduction which 
depends on the magnetic field geometry, is still unknown. 
We have accounted this possibility by allowing the value of 
saturated constant, <f> 3 , to vary in our solutions. Magnetic 
field also has an important role for transferring angular mo- 



mentum along the disks. So the dynamical structure of the 
disk will be affected by magnetic field strength and configu- 
ration. So investigating the magnetized accretion flow with 
thermal conduction is an important issue. 



2 THE BASIC EQUATIONS 

We describe the 2-D hot accretion flow with similar maner 
of Narayan & Yi (1995). We adopt spherical polar coor- 
dinates (r, 0, 4>) for axi-symmetric and steady state flows 

= §i = 0).The fundamental MHD governing equations 
can be written as: 

The equation of continuity gives 







(1) 



The equation of motion 



P^ = -VP-pV0- M V 2 u+( M6 + -U)V(V-u) + -!-JxB(2) 
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The equation of energy 



p[^ + P J^(^)] = Q™ +QB - Qrad + Qcond (3) 



Gauss's law 



V-B = 
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and the induction equation 
V x (u x B) + j?V 2 B 



(4) 
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where p is the density of the gas, p the pressure , e the in- 
ternal energy, u is the flow velocity, B is the the magnetic 
field, J = V x B the current density, r\ the magnetic diffu- 
sivity in which for simplicity it is assumed to be a constant 
parameter (see, e.g., Kaburaki 2000), fi and pb are the shear 
and bulk viscosities . 

The viscous heating rate is defined as the expression: 



Qvi 



2pE ii E ii + {p, b - -aO(V-v) 
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where Eij = ^{vi,j+Vj,i) is a symmetric tensor and is known 
as the rate of the strain tensor. 

We have adopted saturated conduction (Cowie &Mckee 
1977) as: 



Qc 
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where as we have already mentioned F s — 50 s pel is the 
saturated conduction flux on the direction of the tempera- 
ture gradient. Tanka & Menou (2006) have shown that for 
very small <f> s their solutions coincide the standard ADAF 
solutions. They have shown that by adding the saturated 
conduction parameter , <fr s , the effect of thermal conduction 
can be better seen when we will approach to ~ 0.001 — 0.01. 
So, we have investigated the effect of thermal conduction in 
this range. 

So, magnetic reconnection may lead to energy release. 
Also,we can consider the viscous and resistive dissipations 
due to a turbulence cascade. In this study,the resistive dis- 
sipation is defined: 



Q B = -T-J 2 

4ir 



(8) 
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In the right hand side of the energy equation we have: 

Q-\- Q — H~~ Qcond — Qadv — fQ~\- H~~ Qcond 

where Q + = Q vis + Qb, Q- = Qrad and Q adv represents the 
advective transport of energy and is defined as the difference 
between the magneto- viscous heating rate, Q+, and radia- 
tive cooling rate, Q ra d plus the energy transport by conduc- 
tion, Qcond- We employ the parameter / = 1 — to mea- 
sure the hight degree to which accretion flow is advection- 
dominated. When / ~ 1 the radiation can be neglected and 
the accretion flow is advection dominated while in the case 
of small / the disk is in the radiation dominated case. So we 
can rearrange the right hand side of the energy equation to 
fQ+ + Qcond, where / < 1. In general, it varies with r and 
depends on the details of heating and cooling processes. For 
simplicity, it is assumed to be a constant. 

For simplicity,the self-gravity of the disc and the ef- 
fect of general relativity have been neglected. Also, we ne- 
glect radiation pressure in the equations because in optically 
thin ADAFs, P gas > P rad . We adopt the dipolar configu- 
ration for the magnetic field. Also we have neglected the 
(^-component of the flow velocity ug = 0,and the bulk vis- 
cosity of the flow, fit, — 0. Now , we formulate the basic 
equations (l)-(5) in spherical polar coordinates as follows: 



dp 19, 2 x 1 9 , \ 
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dt r 2 dr r < 

The three components of the momentum equations give 
(e.g., Mihalas & Mihalas 1984): r component 

du r u% _ GMp dp 4 dV 8 1 du r 
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The three components of the induction equation 
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Now we have a set of MHD equations which describe the 
dynamical behavior of magnetized ADAFs. The solution of 
these equations are strongly depends on viscosity, resistivity, 
degree of advection and the role of thermal conduction on 
the disks. 

These nine partial differential equations governing the 
non-self gravitating, magnetized advection dominated vis- 
cose flows. These equations relate 15 dependent variables: 
p, p, e, a, p B ,rj and the components of u, J and B. For the set 
of equations, we use the following standard assumptions: 

The kinematic viscosity coefficient, v = ^, is generally 
parameterized using the a-prescription (Shakura-Sunyaev 
1973), 

v = ac s H, (17) 

where H = -Sf- is known as the vertical scale height , 
c s — y/~£ is the isothermal sound speed and the dimension- 
less coefficient a is assumed to be independent of r. So, we 
introduce the parameter r\ as the magnetic diffusivity and 
insert it as a constant parameter in our equations. Both the 
kinematic viscosity coefficient v and the magnetic diffusivity 
r] have the same units and are assumed to be due to turbu- 
lence in the accretion flow. Thus it is physically reasonable 
to express r/ such as v via the a-prescription of Shakura- 
Sunyaev (1973) as follows (Bisnovatyi-Kogan & Ruzmaikin 
1976), 



7/ = IJoCsH. 



(18) 



To determine thermodynamical properties of the flow in the 
energy equation, we require a constitutive relation as a func- 
tion of two state variables. Therefore we choose an equation 



for the internal energy as e = p ^_^ where F is the ratio of 
specific heats of the gas. 

To satisfy V • B = 0, we may introduce a convenient 
functional form for the magnetic field. Owning to the ax- 
isymmetry, the magnetic field can be written as 



B p (r,e) + B <t> (r,e)e 4 



(19) 



Angular momentum is expected to be carried away from 
the disk by magnetic stresses along the externally given 
poloidal magnetic lines of force. In the case of dipole-type 
external field it is transferred to the central accretor (Kabu- 
raki. 2000). The effect of magnetic diffusivity on magneti- 
cally driven mass accretion was studied by Kaburaki (2000) . 
They showed that the effects of resistivity are that magnetic 
field lines do not rotate with the same angular speed as the 
disk matter and thus it suppresses the injection of mag- 
netic helicity and magneto-centrifugal acceleration. So, by 
neglecting the toroidal component of the field, B^, we can 
express the poloidal component, B p , in terms of a magnetic 
flux function ^(r, 6): 

B = B p (r, 6) = -L V x (-^e^ (20) 
27T \rsme } 

It is clear that the basic equations are nonlinear and we 
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cannot solve them analytically. Therefore, it is useful to have 
a simple means to investigate the properties of solutions. 
We seek self-similar solution for the above equations. In the 
next section we will present self-similar solutions of these 
equations. 



3 SELF-SIMILAR SOLUTIONS 

To understand better the physical processes of our viscous- 
resistive ADAF accretion disks, we seek self-similar solutions 
of the above equations. The self-similar method is familiar 
from its wide applications to the full set of MHD equations. 
The self-similar method is not able to describe the global 
behavior of accretion flows, because no boundary condition 
has been taken into account. However, as long as we are not 
interested in the behavior of the flow near the boundaries, 
such solutions are very useful. 

Writing the equations in non-dimensional forms, that 
is, scaling all the physical variables by their typical val- 
ues, bring out the non-dimensional variables. We can simply 
show that the solutions of the following forms, satisfy the 
equations of our model: 



n* 2 = K, where K 



p(r,8) = p a p(e)(r/r )- 3/2 , 
p(r,e)=p P(e)(r/r )- 5/2 , 
u T (r,9) = rfi K (r)l/(0), 
u v {r,0) = r sin 6Q K (r)Q(6), 



(21) 
(22) 
(23) 
(24) 

(25) 
(26) 



where p , p , B a and r provide convenient units with 
which the equations can be written in non-dimensional 
forms. Substituting the above solutions in the equations 
(10)- (16), we obtain a set of coupled ordinary differential 
equations in terms of 9. 
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(29) 
(30) 



Finally, by definition f3 

— o / — ■• 

is an arbitrary constant, so ci = — — 2p R and 

J ' Po V r a ) Po ff 

In equation (30) we have the followings: 

A =plM + mcote ' 



9 ,1 3 , pU 
4 a r) ciP 



9/ 



Equations 27-30 constitute a system of ordinary non- 
linear differential equations for the four self-similar variables 
O, P, Ukp. 

There are many techniques for solving these nonlinear 
equations. Analytical methods can yield solutions for some 
simplified problems. But, in general this approach is too re- 
strictive and we have to use the numerical methods. Here, 
one can employ the method of relaxation to the fluid equa- 
tions (Press et al. 1992) . In this method we replace ordinary 
differential equations by approximate finite-difference equa- 
tions on a grid of points that spans the domain of interest. 
The relaxation method determines the solution by starting 
with a guess and improving it, iteratively. Based on it, this 
system of equations can be solved for all unknowns as a 
function of 9, once we are given a set of boundary condi- 
tions where constraints are placed on the flow. The bound- 
ary conditions are distributed between the equatorial plane, 
= f and the rotation axis, 9 = 0. We can use Narayan 
& Yi (1995)Boundary conditions in both boundaries: the 
boundary conditions at 9 = 



dU dQ. dP dp 







U = 0,p= 0, 



(31) 



d9 d9 d9 d9 

and in this method the boundary conditions at 6 = \ 

are: 

m_dn_dP_dp_ 
d9 ~ d9 ~ d9 ~ d6~ ■ 



(32) 



The boundary conditions on the above equations require 
that variables are assumed to be regular at the endpoints. 
Also the net mass accretion rate (9) provides one boundary 
condition for p: 



p(6)U (6) sin 9d9 = -- 



We obtain numerical solutions for the flows with fixed 
values of r? = 0.1, V = f , a = 0.01, 0.05, 0.1 ,/ = 0.1, 0.3, 0.7 
and <f> s = 0.001,0.007,0.01. We consider ci = 0.8, c 2 = 
2 x 10 3 and (3 = 0.01 (Ghanbari , Salehi & Abbassi 2007) 



4 RESULTS 

We have obtained numerical solutions of equations (27)-(30) 
for a variety of viscosity, a , advection, / and the thermal 
conduction <j> s parameters. The three panels in Figure 1 show 
the variations of various dynamical quantities in terms of po- 
lar angle 9 for with a fixed values of viscosity and thermal 
conduction parameters with a sequence of increasing advec- 
tion parameter /. The top panel displays the dimensionless 
radial velocity U(9). U(9) is zero at 9 = (this is a boundary 
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Figure 1. The self-similar solutions of radial velocity U(8) (top), 
density p(0) (middle) and angular velocity Q(9) (bottom )as a 
function of polar angel Q corresponding to V = 4/3, r]o = 0.1, 
Po = 0.01 and / = 0.1, 0.3, 0.7 for a=0.1 and <j> B = 0.001 



Figure 2. The self-similar solutions of radial velocity U(0)(top) 
and density p{6) (middle) and angular velocity Q(9) (bottom) as 
a function of polar angel 8 corresponding to T = 4/3, rjo = 0.1, 
Po = 0.01 and a = 0.01,0.03,0.05 for / = 0.7 <j> s = 0.001 



condition) and maximum at 9 = ^ . Thus, the inflow ve- 
locity reaches its maximum in the equatorial plane 
and vanishes along the polar axis. As expected , the 
velocity is sub-Keplerian. The middle panel shows the den- 
sity profile p(9) of the solutions. The density contrast in the 
equatorial and polar regions increases with a decrease in the 
advection parameter /. For a given a and 4> s , solutions with 
small values of / behave like standard thin disks, as might 
be expected, these solutions correspond to / — > and so a 
small fraction of energy would be advected. In the opposite, 
/ — » 1, advection-dominated limit, our solutions describe 
nearly spherical flows which rotate far below the Keplerian 
velocity. The bottom panel shows the profile of the angular 
velocity Q(9). Q(9) decreases increase in the advection in the 
accretion disks. We find that in the inner boundary, U(9) is 
essentially independent of advection parameter /. But in in- 
termediate values of 9, the radial velocity is modified by /; 
in the Shadmehri's solutions (2004), two distinct regions in 
the U (9) profile could be recognized. The bulk of accretion 
occurs from equatorial plane at 9 = J to 9 = 9 S , in which 
the radial velocity is zero. While in Narayan & Yi solutions, 
there is no zero inflow in < 9 < -|. Our solutions show 
that in any given 9 the radial velocity is nonzero and when 
we increase the advection parameter the radial velocity will 
be increased. 

Here we may comment about the fact that when the ad- 
vection parameter, /, goes to zero, our disk dose not corre- 
spond to a globally cooling flow, because of the appearance 
of thermal conduction term in the energy transport equa- 



tion. When / goes to 1 our disks are not fully advective, 
because some part of the energy generated by viscosity will 
be transported by thermal conduction. 

Figure 2 displays the behavior of radial and angular 
velocities and density profile for different values of the vis- 
cosity parameter for fixed advection and thermal conduction 
parameters. We find that the value of viscous parameter, a, 
affects quantitatively (but not qualitatively) on the dynami- 
cal variables of the accretion flow. For the larger value of the 
viscous parameter, the radial inflow decreases and the den- 
sity would be increased overally, which is compatible with 
the results presented by Ghanbari et al (2007). 

Three panels in Figure 3 show, for fixed advection and 
viscosity parameters with a sequence of thermal conduction 
parameters,: (1) The radial velocity increases with an in- 
crease in the thermal conduction parameter (2) The density 
profile increases with a decrease in the thermal conduction 
parameter (3) The angular velocity decreases with an in- 
crease in the <j) s in the accretion discs. 

In Fig. 4 we display the isodensity contours in the 
meridional plane. The top, middle and bottom panels dis- 
play the isodensity profiles for different values of conduction, 
advection and viscous parameters, respectively. The panels 
in figure 4 show that the disk seems to be thick. Solution 
with the same / but for different values of a are distin- 
guishable from one another. By adding viscous parameter 
the geometrical shape of the disc becomes thick more and 
more. This advection-dominated solutions have very similar 
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Figure 3. The self-similar solutions of radial velocity U(8)(top) 
and density p(8) (middle) and angular velocity Q(8) (bottom ) as 
a function of polar angel 8 corresponding to T = 4/3, 170 = 0.1, 
Po = 0.01 and <f> 3 = 0.001, 0.007, 0.01 for / = 0.7 a = 0.01 



properties to the approximated solution derived by Narayan 
& Yi (1994), Ghanbari et ql (2007) and Shadmehri (2004). 

The overall structure of the dynamical variables remain 
very close to the original 2D ADAF solutions of Narayan 
& Yi (1995). This solution is denser close to the equator 
than at the pole, and it is some thing like a "thin disk" 
surrounded by a hot coronal atmosphere. 

As the values of tfi s are increased, the solutions substan- 
tially start to deviate from standard solutions, with faster 
radial flow and rotating more slowly (becomes more pressure 
supported) in the equator. By increasing 4> s , density reduces 
in the equatorial plane, and the density profile becomes more 
uniform . It gradually approaches to the spherically symme- 
try. This behavior of the solutions are shared by original 
ADAFs solutions and Ghanbari et al.'s work (2007) 



5 SUMMERY AND CONCLUSION 

The main aim of this investigation was to obtain an 
axi-symmetric self-similar advection-dominated solution for 
viscous-resistive accretion flow with the poloidal magnetic 
field in the presence of thermal conduction. Using the basic 
equations of fluid dynamics in spherical polar coordinates 
(r,8,(p), we have found the self-similar solutions for thick 
discs to derive a set of coupled differential equations that 
govern the dynamics of the system. We have then solved 
the equations using the relaxation method by considering 
boundary conditions and using the a-prescription (Shakura 



Figure 4. isodensity contours for T = 4/3, rjo = 0.1, /3o = 0.01 
and 4> s = 0.001,0.007,0.01 for / = 0.7 a = O.Ol(top) and 
/ = 0.1,0.3,0.7 for a=0.01 and ip s = O.OOl(middle) and a = 
0.01,0.05,0.1 for / = 0.7 <f> 3 = 0.001 (bottom) 

& Sunyaev, 1973) in order to extract some of the similarity 
functions in terms of the polar angle 9. 

We showed that the radial and rotational velocities are 
well below the Keplerian velocity. The Bulk of accretion with 
nearly constant velocity occurs in the regions which extend 
from equatorial plane to a given 9 which highly depends on 
advection parameter /. In a non-advective regime, low /, 
we have a standard thin accretion disk, but for a high / the 
accretion is nearly spherical. 

It is difficult to evaluate the precise picture of radiative 
inefficient accretion flow in the presence of thermal conduc- 
tion with a self-similar method. But this method can repro- 
duce overall dynamical structure of the disks with a set of 
given physical parameters. Even though, conduction heats 
up the accretion flows locally, the reduced density resulting 
from the larger inflow velocity, leads a net decrease in the 
expected level of free-free emissions. The very steep depen- 
dence of synchrotron emissions on the electron temperature 
(e.g. Mahadevan & Quataret 1997) , suggests that hotter so- 
lutions (with conduction) maybe more efficient radiatively 
( Tanaka & Menou 2006). From self-similar solutions alone, 
we can determine how the global structure of the flow can 
be affected by thermal conduction. For small enough val- 
ues of 4> s the solutions remain very close to the Ghanbari 
et al. (2007) solutions. The main difference between their 
solution with standard solution of 2D ADAFs of Narayan & 
Yi (1995a) is the presence of dipolar magnetic field and its 
correspond resistivity. In our case, we add extra physics by 
adding the thermal conduction as a mechanism for energy 
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transport. By adding the thermal conduction the solution 
start deviating substantially from the original ADAFs, with 
faster radial inflow at the equator. These results well agree 
with Tanaka & Menou (2006). 

The presence of magnetic field with the poloidal con- 
figuration will affect the role of thermal conduction. Com- 
pare to non-magnetic field solution, (Tanaka & Menou 2006) 
in our case the existence of magnetic resistivity can pro- 
duce more energy to be advected. The B-field configuration 
can also affect the energy transportation along the accretion 
disks. However, the main aim of this work is to study the 
quasi-spherical magnetized flow, directly by solving the rel- 
evant MHD equations. Although, we have made some sim- 
plifications in order to treat the problem analytically, our 
self-similar solutions show that the input parameters, such 
as thermal conduction and viscose parameters, magnetic 
field and its resistivity can really change typical behavior 
of the physical quantities of the ADAF disks. Of course, our 
self-similar solutions are too simple to make any compar- 
ison with observations. But, we think one may relax self- 
similarity assumptions and solve the equations of the model 
numerically. This kind of similarity solutions can greatly fa- 
cilitate testing and interpreting the results. 
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